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ABSTRACT 


The effects of thermal imbalance (S-^ ^ 0) on stellar pul- 
sational stability have largely been ignored in the literature. 

A preliminary investigation of such effects is made here using 
the linear, quasi-adiabatic pulsation theory. Analyses of the 
standard model and of a white dwarf show that, for these cases, 
"ordinary" terms in the stability integrals greatly outweigh 
the "extra" terms considered here. The latter become more 
important when substantial ionization zones exist in the stellar 
matter. It is argued that the influence of thermal imbalance 
on pulsational stability should be quite small for stars cross- 
ing the H-R diagram in early post-main sequence evolution, and 
for cooling degenerate stars. On the other hand, during pre- 
main sequence contraction the thermal imbalance terms are more 
likely to be important, while for thermal runaways in shell- 
burning stars, these terms are almost certainly crucial. 
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I • INTRODUCTION 


In the familiar linear quasi-adiabatic theory of pulsations, 
developed in general form by Thomas (1930) , the rate of change 
of the dynamical energy of oscillation E may be written (Ledoux 1958) 
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where the bar indicates a time average over the pulsation period , 
and the subscript 1 denotes equilibrium quantities. U and S are, 
respectively, the specific internal energy and entropy of the gas, 
p is the mass density, and the quantity 6 preceding a variable 
indicates a first-order departure from equilibrium. 


In the stellar pulsations most commonly studied the equilibrium 
entropy change 

• 1 SL r 

S, = - [e- ] / 

1 T am 


where e is the nuclear energy generation rate and L r the luminosity 
at distance r from the stellar center, is negligible and the right 
side of equation (1) reduces to the first integral with 
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For those epochs of stellar evolution during which S^O 
(thermal imbalance) , the expression (2) must be rewritten as 
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and contributions from the second, and perhaps the third, in- 
tegral on the right side of (1) considered. The latter integral 
was ignored by Thomas as vanishing in the mean, but was included 
by Ledoux who pointed out that second order contributions from 
it might subsist. 

The purpose of the present work is to consider both the second 
integral on the right of (1), and the extra term introduced in the 
first integral due to the replacement of (2) by (3) . Because adia- 
batic evolution (Sj =0) is so often a good approximation in stellar 
model calculations, these terms have, with few exceptions, been 
ignored in the literature since the general review of Ledoux (1958) . 
In what follows, we shall attempt a preliminary investigation of 
the problem, proceeding through the use of simple stellar models 
and semiquantitative estimates to get an idea of the relative size 
of the thermal imbalance terms and to draw some conclusions regard- 
ing their importance in various phases of stellar evolution. 


In section II we write equation (1) in more convenient form. 
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Sections III and IV are devoted to evaluating the stability of 
the standard model and that of a cooling white dwarf, respec- 
tively. In section V we consider what effects the presence of 
ionization zones in the stellar matter might have, and in section 
VI discuss the results of the various calculations. The final 
section considers two cases in which thermal imbalance terms are 
likely to be important, and for which further investigation should 
prove fruitful. 


II. THE STABILITY INTEGRAL 


Equation (1) gives the rate of gain or loss of pulsational 
energy over a cycle. Let us define the first two integrals on 
the right-hand side as 
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This represents the contribution to the pulsation analysis of the 
terms we will consider. 

In the usual notation we have 


-^2 = 2 ( r ,- i ) , 
3Sdp p ^ 



and using (3) , we may rewrite (4) as 
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Let us now make the assumption that no sources or sinks of 
subatomic energy exist in the stellar matter. In that event, 
defining the pulsational quantities 



6£. r 
p 



and with the adiabatic condition 


t = (r 3 -l)2' 


we have 



(r.-l)^ i 6L ) dq 
3 L dq r 


(6) 


+ f 1 (r -l) 2 -M z 2 dq - fjd 3 3 u m g 2 

0 3 3q a 2T , sap 2 ^ 


where f = L r /L and q = m/M. 


For future convenience we shall rewrite (6) as 
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L 


=I 1 +I 2 + J 2' 


where the integrals (including signs) are labelled in order. 

We shall call terms coming from the first integral in (6) (subscript 1) 
the "ordinary" terms, while those which arise from the other two 
integrals (subscript 2) will be referred to as "extra" or "thermal 
imbalance" terms. 

A star is considered pulsationally stable if L'/L <0. In the 
event that contributions from subatomic energy do exist in the star, 
one must evidently return to equation (5) . 


III. THE STANDARD MODEL 

Following Eddington (1959) , let us define the quantity ^ by 
the relation 

f “ rjq • 

Then, if k represents the opacity, the standard model is obtained 
for 

■qH = const. 

For our analysis, we shall choose the combination 

•q = const, k = const. 

In that case, since f = q = 1 at the stellar surface, we must have 
r\ = 1 , and our choice corresponds to a physically reasonable case 
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of electron scattering opacity, and luminosity increasing linearly 
outward with mass fraction. 

The standard model is a polytrope of index 3 with a constant 

ratio p of gas pressure to total pressure through the star. The 

2 

value of p depends only on the product Mq , where q is the mean 
molecular weight of the matter, and may thus be chosen arbitrarily. 
All of the non-pulsational quantities in (6) may be evaluated with 
the aid of Emden tables (British Association for the Advancement of 
Science 1932) . 

The relative pulsational amplitudes x, z of the standard model 
have been given by Schwarzschild (1941) for various values of p, 
both in the fundamental and higher modes. We have chosen the lowest 
of these values of p in the fundamental mode as the combination most 
conducive to pulsational instability (see, e.g., Simon and Stothers 
1969) . This is 

p = 0.510, r 3 - 1 = 0.372. 

Furthermore, we may write 


1 

L Sq 


(6L r ) 


M + % — 
3q 3q 


where, for a radiative star with constant opacity, 


l = 4(x + t) 
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dt 
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The remaining unexplicit quantity in (6) is 
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we easily find 


= (1 - 4(r, - i) L 

dp's p • 1 3 
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and after some algebra 


2 
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2t asap 2 


w x (p) = 


327p -1128f7+ 1312 p- 512 

9(8-7p) 3 


The function W-^ ( p) is always negative, with endpoint values 
of -0.111, and a minimum of -0.149 at p = 0.93. A glance at the 
last integral in (6) shows that so long as the luminosity increases 
(decreases) outward, this term serves to energize (damp) pulsa- 
tions. For p = .510, we have W^ ( p) = -0.119. 

We are now in a position to evaluate Lp/L . Results are dis- 
played in Table 1. It is easily seen that contributions arising 
from the additional term in 6S (I 2 ) and from the integral J 2 are 
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approximately equal, and both contribute to energizing the pulsa- 
tion. Their combined effect, however, is quite small, amounting 
to less than 8% of the ordinary damping. 


IV. WHITE DWARF 

We shall discuss here the model of Marshak (1940) treated by 
Ledoux and Sauvenier-Gof f in (1950) . Though somewhat outdated, 
this model is quite adequate for our purposes. In particular, 
the luminosity distribution f(q) is rather similar to that obtained 
from more modern treatments of cooling white dwarfs (e.g. Vila 
1969) . 

Following Ledoux and Sauvenier-Gof fin , we shall take 

z = -3x = const. 

t = (r T -i) z 

where 

_ x 2 + 2 

r T " 1 ~ 3(x^+l) 

and x is the ratio of the Fermi momentum and me. 

We begin with the second integral in (6) . Integrating by parts 

o 2 1 

I 2 = 4x - 18x J* (r T -i) fd(r T -l) , 


we obtain 
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where we have taken r -1 = 2/3 at the stellar surface. Given the 

rri 


model parameters we may easily evaluate the integral, obtaining 

X 1 C ~ 2 

I = 1 .6 x . 

2 

Turning now to we must calculate 
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Using the familiar expressions for a degenerate electron gas 
(Chandrasekhar 1939) and retaining only terms of lowest order in 
(kT/mc 2 ) , we find 


dx | _ x 
dp B 3p 


and finally, 

2 ^3 rT 

£ i-U s W,(x) 

2T dSdp 2 ^ 


2x 4 + 7 x 2 + 2 . 
18 (x 2 +l) 2 


This function falls to a minimum of -0.153 at x=l . At the endpoints 
(x=0, x — >co) , W 2 (x)= -0.111. Since W 2 (x) < 0, it turns out again 
that, for luminosity increasing outward, is an energizing term. 

We have for the present case 

1 

J = -9x 2 J W (x)df. 

2 0 
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The values of x for our model range from 2.4 at the center to 0.4 
at the surface, encompassing values of j (x) | in the narrow range 
0.13-0.15. Using an average value |w 2 (x) |= 0.14, we obtain 

~2 

J 2 =1.3 x . 

The remaining term in (6) is given by Ledoux and Sauvenier- 
Goff in: 1^ = -18x . 

Setting x = 1, we again summarize our results in Table 1. 

In the present case, as with the standard model, the contributions 

= 1.6 and J 2 = 1.3 both serve to energize, and are roughly equal. 

Here, however, their effect is somewhat more important, offsetting 

about 16% of the ordinary damping. 

V. IONIZATION ZONES 

When ionization zones cover substantial regions of the stellar 

matter, the pulsational stability of a star can be strongly affected. 

In the absence of thermal imbalance, such effects have been studied 

in the literature in great detail. 

When thermal imbalance exists in a star, the terms I~ and j 

z* z. 

begin to contribute, and they are affected by ionization through 
the thermodynamic coeffients (r 3 ~ 1) and ( p 2 /2T) S 3 U/SS dp 2 . (It should 

be noted that the quasi-adiabatic theory may not be adequate here, 
depending upon the location of the ionization zones in a given star. 
The theory should be good enough, however, to give us some picture 

of the relative importance of thermal imbalance to the pulsational 
stability) . 
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If in the presence of ionization zones, we impose the restriction 
that only one critical state of ionization exists in a given range of 
temperature and density, then we may write (Ledoux 1958) 


[16-12 p— 1 • 5 p 2 + p (4—1 • 5 p+ ph) Aip p J f 1 +A p ip J 
r l ~ 12-10 .5p + p(l .5+h) A t 
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Further, if X represents the abundance (by number) of the 
element with critical ionization, y the fraction of the electrons 
in question that have been removed, and x the total number of free 
electrons per ion, then 

B = Xx^llz- Z) 

x(l+x) + Xy(l-y) 


The quantity h is the ratio of the ionization potential to kT: h=.r/kT. 
Once more, we must calculate the quantity 


2T SSBp 


7 = w 3 (p,T) = %(r 3 -i) [(r 3 -i)-i] + % 


S-inp 


S (r 3 _1) 


( 8 ) 
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The necessary derivatives are (dp/dp)g, given once more by 
(7) , and 


5x | _ (1+x) 

dp S p 


C A p ,t + tr 3 -i) a t , p 1 • 


Finally, after cumbersome but straight-forward computation we obtain 
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To get some idea of the size of W 3 (p,T) we have calculated it 
for a range of densities and temperatures for four different cases 
of critical ionization: HI, He I , He II, C VI . The ratio of parti- 
tion functions for the states involved was always taken to be the 
ratio of statistical weights. The latter quantities, along with 
values for ionization potentials, were taken from Unsold (1955) . In 
each case a composition consisting solely of the element in question 
was assumed for simplicity. This means, of course, that the effect 
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of the ionization zones will be overestimated. 

Table 2 presents these results. The last entry in the table 
is the ratio 

W 3 (p,T) 

R = 

( r 3-i) 2 

Since ordinary damping terms (at least for the case of radiative 

2 

damping) as well as the integral I 3 will be proportional to (r^-l) , 

the quantity R gives some guide as to the importance of the last 
integral in (6) . 

In evaluating Table 2, it will be useful to compare W 3 (p,T) 
with the quantities W ( p) and W (x) calculated in previous sections. 

J- ^ 

We first note that W^ is not necessarily negative. In certain 
regions, the derivative [a(r^-l) /3-tnp ] can become positive and 
large enough to outweigh the first term in (8) . This happens in 
general for large positive Q and for 3-1. In that case, a star 
with luminosity increasing (decreasing) outward would tend to be 
damped (energized) by the last integral in (6) . However, this effect 
seems to occur in relatively few parts of the (p,T) plane. Overall, 
those regions with W^ < 0 are likely to dominate, giving a result 
qualitatively the same as in the absence of ionization zones. 

Perhaps a more important difference involves the size of the 
effect. We see from Table 2 that W^ can become quite large in ab- 
solute value, exceeding the maximum values of J W-^ | and | W 2 | by a 
factor as great as 3. Furthermore, the ratio R attains large negative 
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values (e.g., He I, log p = -8, log T = 4.2) , particularly in 
low density regions. The largest values occur for p -1 < 0.1, and 
run as high as R =-5.78 for the cases calculated. (As one goes to 
ionizations with larger potentials, e.g. C VI , the effect begins 
to diminish. This is due to increasing domination by radiation 
pressure in the low density regions, and to the fact that the ratio 
x of electrons to ions is increasing. In fact, as g -> 0 or x -> 

W -> W ) . The corresponding ratios for W-, and W can never exceed 

•J X ■*“ ^ 

unity in absolute value, and are generally much smaller. 

Thus it turns out that the thermodynamics of ionization zones 
is such that the contribution to stability analysis of "thermal 
imbalance" terms may become greatly enhanced. Whether this actually 
happens depends of course on the detailed properties of the stars 
in question. We shall have more to say about this in the next section. 

VI. DISCUSSION 

Up to now we have seen that the thermal imbalance terms investi- 
gated will, in general, tend to energize pulsations in stars with 
luminosity increasing outward. This will be the case for gravita- 
tionally contracting stars or stars experiencing thermal cooling. 
Expansion, on the other hand, will tend to damp pulsations. 

In sections III and IV it has been shown for two different cases 
that the extra terms make a small contribution compared to that from 
ordinary damping terms. It is simple and instructive to further 
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compare these extra terms with the energizing due to nuclear re- 
actions. For the case of the white dwarf, Ledoux and Sauvenier-Gof f in 
(1950) have shown that nuclear energizing would be enough to overcome 
the damping for very modest values of the temperature exponent v 
(v = 9.5 for a core source; \, = 2.6 for a shell source) . Thus, this 
energizing is at least six times as effective as the thermal im- 
balance contribution given in Table 1. 

For the standard model, a comparison is not so straight-forward. 
However, since, apart from the temperature exponent, the amount of 
nuclear energizing will depend mainly on the size of pulsational 
amplitudes in the burning region (let us say in the core) , it seems 
reasonable to compare the standard model to a model with a similar 
value of relative radius amplitude at its center. One such model 
in the literature is that of a 28.2 Mq main-sequence star ( Schwarzschild 
and Harm 1958) . It has a central radius amplitude x c = 0.39, while 
the value for standard model of section III is x c = 0.32. With a 
temperature exponent v = 13 , Schwarzschild and Harm (1959) found nuclear 
energizing L pN /L = 4.4 — a value more than an order of magnitude 
larger than that due to thermal imbalance in the standard model 
(Table 1) . 

The above comparisions are, of course, somewhat artificial 
since 1) any region of a star in nuclear thermal equilibrium (S^~ 0) 
cannot contribute to the thermal imbalance terms, and 2) the standard 
model could not represent in detail a star with physically reasonable 
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nuclear processes. Nevertheless, we may safely say that, given the 
thermodynamics of sections III and IV, ordinary pulsation terms will 
tend to greatly outweigh thermal imbalance terms, with the latter 
providing only relatively small corrections. 

For the thermodynamics of section V (ionization zones) the 
outlook changes somewhat. We have seen in this case that the coeffi- 
cient jw 3 | can become relatively large, tending to enhance the con- 
tribution due to thermal imbalance. On the other hand, such effects 
occur mainly in diffuse regions, where the density is low. The 
classic example is in the ordinary Cepheid variables. In the well- 
known modern calculations of Cepheid variability (Christy 1966ei, 

Baker and Kippenhahn 1965) , thermal imbalance is ignored, and the 
equilibrium luminosity taken as constant in the energizing region. 

Due to the extreme central condensation of Cepheid models, only the 
outer 1 or 2 per cent of the stellar mass is affected by the pulsa- 
tion (Christy 1966b) , and even if thermal imbalance does exist in 
these layers, it is hardly to be expected that they contribute 
enough luminosity to sensibly affect the stability. 

Further, it seems possible on the basis of our results to 
rule out important contributions for any star expanding or con- 
tracting across the H-R diagram during normal post-main sequence 
evolution, at least through core helium burning. Although con- 
siderable regions of such stars can be in a state of thermal im- 
balance (see, e.g., Iben 1966) , the central condensation of these 
objects (P c /<p> > 10 3 - 10 5 , as compared with p c /< p > = 54 for the 
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standard model) is such that radiative damping will almost certainly 
crush any of the extra terms we have considered. 

Similarly, cooling degenerate bodies such as white dwarfs 
are highly unlikely to become pulsationally unstable due to contribu- 
tions from thermal imbalance, such contributions being simply too 
small to overcome damping in the thin radiative zones at the surface 
of these stars. 


VII. AREAS FOR FUTURE WORK 

We begin this section by noting that our investigation remains 
incomplete due to the omission of the last integral in (1) . However, 
there are reasons to believe that inclusion of this term will not 
change qualitatively the conclusions of the previous section. Call- 
ing the integral in question LK 2 , taking the time average and nor- 
malizing with the luminosity, we obtain 


K 


2 



where (6p/p) 2 indicates the surviving time-averaged second order 
amplitudes . 

We have already seen in the cases studied that the integrals 

I 2 and J 2 provide small, approximately equal contributions. Since 

. ~2 

we must expect that (6p/p) 2 « Z , it follows that K 2 » I 2 , J 2 , or 
perhaps somewhat larger due to the relative largeness of (p 0 - 1) 
compared with the thermodynamic coefficients of I 2 and J 2 . Even 
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if the integral K 2 proved to have the same sign as I 2 and J ^ in 

the cases investigated, it would need to be nearly ten times larger 
for the white dwarf and more than twenty times larger than I 2 or 

J 2 for the standard model in order to qualitatively affect the 

stability. Although a full second order theory must be developed 

to exactly evaluate the size of K 2 , such large values seem unlikely. 

On the other hand, when the terms I 2 and J 2 begin to become 
comparable to the ordinary damping or energizing terms, the in- 
tegral must obviously be considered in detail. We shall suggest 
two stages of stellar evolution for which this may be the case: 

1 . Thermally unstable shell-burning stars . Rose has found 
that certain thermally unstable models become pulsat ionally unstable 
as well, for both helium (Rose 1967) and hydrogen (Rose 1968) shell 
burning. The seat of the instability is strongly enhanced nuclear 
energizing arising from a thermal runaway in the shell. However, the 
flood of photons released in the shell is almost totally absorbed in 
the expanding layers above, with the luminosity dropping by orders 
of magnitude from shell to surface. Although Rose has not published 

the runs of luminosity or pulsational amplitudes, he does give 
enough information to enable us to put limits on and for his 
unstable models. 

Table 3 shows the relevant quantities for two typical models: 

3B (Rose 1967) and 2A (Rose 1968) . The entries are, in order, the 
shell and surface luminosities, the relative pressure amplitudes 
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p = 6P/P, the nuclear energizing (normalized by L) , and the 

ordinary damping 1^ . The last two entries give limits on the 

thermal imbalance terms calculated as follows. 

Since the envelopes of Rose's models are non-degenerate or 

2 

only slightly degenerate, we have used (r\- 1) and W ( p) as the 

3 X 

thermodynamic coefficients for and , respectively, and have 

set p = 1 for simplicity. Using these quantities and neglecting 

f ,, with respect to f , n n , we obtain 
surface r shell 


+ h ■ 


1 f 

5 shell 


< P > 


(9) 


where < p > is an average over the envelope luminosity distribu- 
tion, and = L„A- 

The sum (9) is negative (i.e., a damping term) as we expect 

for luminosity decreasing outward. Using the limiting values of 

p as given by Rose, we obtain the last two entries in Table 3. 

Consider the minimum damping due to thermal imbalance (I„ + J ) 

2 2 min 

For model 2A it is comparable to the ordinary damping 1^ , while for 

3B it far exceeds I . In both cases the nuclear energizing dominates, 

so there is no qualitative effect. On the other hand, if we use 

the upper limit (I„ + J ) , the thermal imbalance damping completely 

z 2 max 

overwhelms all other terms and the instability is nullified. Because 
pulsational amplitudes will tend to drop off rapidly from the sur- 
face inward, the true value of I ^ + J ^ is probably closer to the 
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lower than to the upper limit. However, the size of even the minimum 
terms indicate that thermal imbalance contributions must be taken 
into account in any evaluation of the pulsational stability during 
a thermal runaway. To do this will require a full second order 
theory . 

2. Stars in pre-main-sequence contraction . Such stars are 
favored by three characteristics: 1) a large percentage of the 

luminosity is provided by gravitational contraction; 2) the central 
condensation is lowered substantially by the influence of convec- 
tion; and 3) deep ionization zones can exist in such stars, depend- 
ing on mass and the state of evolution. Detailed pre-main- sequence 
contraction models have been constructed by Iben (1965) and Ezer 
and Cameron (1967) . 

A start on a stability analysis for such stars has been made 
by Kato and Unno (1967) and Okamoto (1967) . The former authors 
developed a second order theory and , with the aid of a number of 
restrictive assumptions, managed to write the stability integrals 
(1) in approximate form. This form was in turn used by Okamoto to 
test the stability of a fully convective homologously contracting 
polytrope of index 1.5. Ionization zones were neglected. 

Emphasizing the tentative nature of the calculation, Okamoto 
concluded that stars with masses < 2M Q were probably pulsationally 
unstable against energization due to thermal imbalance. Establish- 
ment of this conclusion must await the development of a full second- 
order theory including all the terms in equation (1) - 
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TABLE 1 


Values of the Stability Integrals 


Standard Model 


*1 

I 

2 

J 2 

Lp/L 

-4.8 

0.20 

0.18 

-4.4 


White Dwarf 


-18 


1.3 


-15 















TABLE 2 


Some Thermodynamic Properties of Ionization Zones 


H I , *(ev) = 13.59 


log p = -9.0 


log p = -5.0 


log p = -2.0 


log T 

3.7 

3.8 1 

4.0 

m 

■ 

KB 

mi 

— 

4.1 

4.4 

4.5 

f 

4.6 

4.2 

4.6 

5.3 

y 

0.000 

0.005 



0.523 

0.956 

0.997 

0.001 


0.045 

0.792 

0.943 

0-984 

0.006 

C .219 

0.945 

8 

0.996 


0.980 

0.970 

0.943 

1.000 

QQ9 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1-000 

r i - 

1.625 





1.602 

1.430 

1.277 

1.339 

1.499 

1.620 

1.560 

1.453 

1.660 

r 3 - 1 

0.620 

0.319 

wm 


0.525 

0.597 

0 .407 

0.231 

0.261 

0.449 

0 .602 

0.544 

— 

0-372 

0.648 


-0.300 

-0.462 

+0.008 

+0.287 

+0.150 

-0.290 

-0.322 

-0.060 

+0.108 

+0.304 

+0.180 

-0.149 

+0 .058 

+0 .022 

w (p,T) 

-0.268 

-0.340 






-0.119 

-0.042 

+0.028 

-0.030 

-0.199 

-0.088 

-0.103 

1 1 

-0.70 

-3.33 

-4.80 

+0.92 i 

-0.18 



-0 .74 

i 

H* 

*~J 

O 


-0.62 

+0.14 

-0.08 

-0 .67 

-0.64 

-0.24 

| 

He I, x{ev) - 24.58 | 

log p = ~8*0 

log p - -4.0 

log p = -1.0 


log T 

4.0 

4.1 


m 

IB 

4.. 2 


4.7 

4.8 

4.4 

4.6 

1 

4.8 

5.0 

5.5 

y 

0 .002 

0.036 




0 .004 

0.479 

0.962 

0.991 

0.005 

0.061 

0.283 

0.646 

0.979 

8 

0.988 

0.977 


ESI 


1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1 .000 

h 

1.465 

1.150 

1.099 

1.192 

1.527 

1.463 

1.230 

1.508 

1.631 

1.547 

1.401 

1.407 

1.505 

1.663 

r 3 - 1 

0.451 

0.129 

0.077 

0.165 

0.463 

0.446 

0.157 

0.466 

0.619 

0.531 

0.346 

0.321 

o 

to 

o 

0.657 

B 1 (r,-D 

-0.692 

-0.054 





-0.025 

-0.416 

+0.026 

+0.367 

+0.161 

-0.194 

-0.022 

+0.056 

+0.112 

+0 .013 

W„(p,T) 

■ j 

-0.470 


-0.034 

-0.006 

-0.137 

-0.331 

-0.054 

+0.059 

-0.038 

-0.222 

-0.124 

. 

-0.081 

-0 .066 

-0.106 

R 

-2.31 

-5.00 

-5.78 

-0.22 

-0.64 

-1.66 

-2.16 

+0.27 

-0.10 

-0.79 

-1.04 

-0.78 

-0.37 

-0.25 































































log p = -8.0 



1 

■ 

1 

log T j 

4.4 

4.5 

4.6 

4.7 

Y 

0.000 

0.019 

0.524 

0.969 

0 

0.912 

0.840 

0.766 


r i 

1.565 

1.237 

1.041 

1.353 

V 1 

0 .516 

0.247 

0.120 

0.331 


-0.102 

-0.366 


m 

W 3 (p,T) 

-0.176 

-0 .276 

-0.051 

-0.003 

R » 

j 

-0.66 

-4.54 

-3.53 

-0.03 




log p = “7 .0 

log T 

5.4 

5.5 

5.6 


y 

0.002 

0.206 

0.928 

0.997 

i 

P 

0.094 

0.051 

0.029 

0.015 

r i 

1.348 

1.304 

1.333 

1.336 

v 1 

0.337 

0.324 

0.333 


ipls ( V« 

-0.003 



0.000 

w 3 (p,T) 

-0.113 

-0.122 

-0.108 

-0.111 


R 


-1.00 


-1.17 


-0.97 


- 1.00 


TABLE 2 (Continued) 
















































































































































TABLE 3 


Pulsational Characteristics for Two Thermally Unstable Models 
3B (Rose 1967) and 2A (Rose 1968) . 



3B 

2A 

l n / l © 

1.5 X 10 6 

o 
1 — 1 

X 

CM 

L/L© 

110 

200 

p ( surface) 

-12 

-15 

p (shell) 

-0 .29 

-0 .50 

N 1 

3.5 x 10 3 

50 

*1 

-4.5 

-6.3 

(I 2 + J 2> min 

-2.4 x 10 2 

-5.0 

(I 2 + J 2 >max 

-4.0 x 10 3 

-4.5 x 10 3 


